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Abstract 

We present an algorithm that determines the Galois group of linear difference equa¬ 
tions with rational function coefficients. 


1 Introduction 

The current algorithms for computing the Galois group of linear difference equations were 
only valid for the equations of special types, such as the second order equations, the equa¬ 
tions of diagonal form or with constant coefficients and so on. In [9], a difference analogue of 
Kovacic’s algorithm was developed for linear difference equations of order two. In [23], algo¬ 
rithms for linear difference equations of diagonal form were developed. For linear difference 
equations with constant coefficients, an algorithm can be found in [22], where the author 
further showed that there is a recursive procedure that derives the Galois group from the 
ideal of algebraic relations among solutions, and vice versa. In [15], Maier gave upper and 
lower bounds for the Galois groups of Frobenius difference equations over (F q (s, t), 4> q ), where 
4> q (s ) = s q and cf> q (a ) = a for all a € F q (t). On the contrary, algorithms for computing the 
Galois groups of linear differential equations have been well-developed (see El tia uni nu). 
Particularly, in HU. Hrushovski developed an algorithm that calculates the Galois groups 
of all linear differential equations with rational function coefficients. His algorithm involved 
many arguments from logical language and has recently been reworked by Rettstadt in m 
and by the author in [7J. Here, in this paper, we develop an algorithm for computing the 
Galois group of linear difference equations with rational function coefficients of arbitrary 
order. Our algorithm can be considered as a difference analogue of Hrushovski’s algorithm. 

The philosophy of computing the Galois groups of linear difference equations is quite 
similar to that of linear differential equations. The Galois groups of these two kinds of equa¬ 
tions are linear algebraic groups over the field of constants. Hence bounds for the defining 
equations of linear algebraic groups developed for the differential case can be applied to the 
difference case without any modification. However, there exist some results in differential 
algebra whose difference analogues are not correct any more, and vice versa. For example, 
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associated primes of a radical differential ideal are again differential ideals, while those of 
a radical cr-ideal need not be cr-ideals but o^-ideals for some integer 6. This forces us to 
consider cr 5 -ideals. Another example is that the Picard-Vessiot extension ring for linear dif¬ 
ferential equations is not necessarily the coordinate ring of a trivial torsor for the Galois 
group, while that for linear difference equations is the coordinate ring of a trivial torsor. 
This implies that one only needs to consider objects such as hyper geometric elements that 
are defined over the basic field. 

Throughout this paper, k stands for the field of rational functions in x with coefficients 
in Q, the algebraic closure of the field of rational numbers, and k stands for its algebraic 
closure. The difference field which we are interested in is the field k with an automorphism 
a given by cr(x) = x + 1 and <r(c) = c for c € Q. Consider the following linear difference 
equations 

a(Y) = AY (1) 

where Y is an n x 1 vector with indeterminate entries and A € GL n (/c). Let R be the Picard- 
Vessiot extension ring of k for P). The Galois group of (JTj) over k, denoted by Gal (R/k), is 
defined to be the set of u-fe-automorphisms of R, i.e. ^-automorphisms of R that commute 
with a. Let F be a fundamental matrix of P with entries in R, i.e. F £ GL„(i?) satisfying 
(j(F) = AF. Then for any </> € Gal (R/k), 4>(F) is another fundamental matrix of (P). Thus 
there exists [(f)] € GL n (Q) such that (f>(F) = F\<p]. The map given by 0 —>• [(f)] is a group 
homomorphism of Gal(R/k) into GL n (Q). Denote by G the set {[0] | cf) £ Gal(R/k)}. It 
was proved in (Theorem 1.13, page 11 of [23]) that G is a linear algebraic group defined over 
Q. The reader is referred to Chapter 1 of [23] for more information about the Galois theory 
of linear difference equations. 

The group G can be reformulated as the stabilizer of some ideal in a cr-ring, which we 
describe below. Let Y denote an n x n matrix (ytj), where the ijij are indeterminates. 
Sometimes, in brief, we also consider Y as a set of indeterminates. By setting cr (Y) = AY, 
one can extend a from k to k[Y, l/det(Y)] so that it becomes a difference extension ring of 
k. The results in Section 1.1 of [23] imply that R is isomorphic to k\Y, 1/det(Y)]// for some 
maximal cr-ideal I. Define an action of GL n (Q) on k[Y, 1/det(Y)] given by g ■ Y = Yg for 
all g £ GL n (Q). Suppose that J is an ideal of k[Y, 1/det(Y)]. The stabilizer of J, denoted 
by stab(J), is defined as 

stab(J) = {g € GL n (Q) | P(Yg) € J, V P € J}, 
which is an algebraic subgroup of GL n (Q). Set 

I F = {P £ k\Y,\/ det(Y)] | P{F) = 0}. 

Then Ip is a maximal cr-ideal and G = stab(Jir). By the uniqueness of the Picard-Vessiot 
extension ring of k for (P), one sees that for any maximal cr-ideal I of k\Y, 1/det(Y)], there is 
g € GL n (Q) such that g-1 = Ip. From this, one can readily verify the stabilizers of maximal 
cr-ideals in k[Y, 1/det(Y)] are conjugated. In other words, as linear algebraic groups, these 
stabilizers are isomorphic. Therefore we shall also call the stabilizer of a maximal cr-ideal 
of k\Y, 1/det(Y)] the Galois group of P) over k. Using the Grobner base method, one can 
obtain the defining equations of stab(J) easily once a Grobner basis of I is known. Therefore, 
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the above definition indicates that finding a maximal cr-ideal of k\Y, l/det(Y)] will suffice 
to determine the Galois group. We shall give in this paper an algorithm that computes a 
maximal cr-ideal of k[Y, 1/det(Y)]. 

The rest of the paper is organized as follows. In Section 2, we introduce some basic 
results that provide the theoretical background of our algorithm. Meanwhile, we introduce 
some basic definitions such as proto-groups, proto-maximal cr-ideals and so on. In Section 
3, we show how to compute a proto-maximal cr-ideal. In Section 4, we describe a method 
to extend a proto-maximal cr-ideal to a maximal cr 5 -ideal so that one can easily obtain a 
maximal cr-ideal by taking the intersection of ideals. In Section 5, the methods developed 
in the previous sections are summarized as an algorithm, and an example is presented to 
illustrate the algorithm. In Appendix A, we describe a method to find coefficient bounds 
for generators of a proto-maximal cr-ideal. In Appendix B, an algorithm for computing cr* 5 - 
hypergeometric elements in k[Y, l/det(Y)]//i rr is developed, where I irr is a prime cx^-ideal. 

2 Some basic results 

In this section, we shall introduce some basic results about proto-groups, A;-torsors and 
several related problems whose algorithmic solutions will be needed in our algorithm. 

2.1 Proto-groups 

As in the differential case, bounds on algebraic subgroups of GL n (Q) play a central role in 
the main algorithm presented in this paper. Let H be an algebraic subgroup of GL n (Q). 
For the ease of notation, we shall use H(k ) (resp. H(k )) to denote fc-points (resp. ^-points) 
of H. We shall say H is bounded by a positive integer d if there is a set § C Q[Y] such that 
H is the set of zeroes of § in GL n (Q) and elements of § are of degree not greater than d. In 
brief, H u stands for the algebraic subgroup of H generated by unipotent elements and H° 
denotes the identity component of H. 

Definition 2.1 Let G,H be two algebraic subgroups o/GL n (Q). H is said to be a proto¬ 
group of G if they satisfy the following condition 

H u < G° < G < H. 

In the case that G is the Galois group of m over k, H is called a proto-Galois group of (QP. 

Remark 2.2 Suppose that H is a proto-group of G and H is an algebraic subgroup satisfying 
G < H < H. Since H u < H u , one sees that H is also a proto-group of G. 

For the convenience, we introduce the following definition. 

Definition 2.3 A cr-ideal I in k\Y, 1/det(Y)] is called proto-maximal if stab(I) is a proto- 
Galois group of W- 

The key point of Hrushovski’s algorithm is the following proposition, which is also the core 
of our algorithm. 
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Proposition 2.4 (Corollary 3.7 of ! Corollary B.15 of w One can find, an integer d 
only depending on n such that for any algebraic subgroup G of GL n (Q), there is a proto¬ 
group of G bounded by d. Particularly, given linear differential equations, there exists a 
proto-Galois group of it bounded by the integer d. 


The integer d can be explicitly given as follows (see Corollary B.15 of [?] for details). Set 


= max < 


m 2 +(2n) 3 ' sn \ 
1 n 2 > 


3 on 2 /n 2 + (2n) 3 ' 8 " 
«2 = «r(2 nf 8 V 

\ n z 


(2) 


and 


, 2 1 \ j f K i + 1\ 1 7Y', t ( \ (n 2 + 1 

n 3 = ^ 2(^1 + 1) max ^ [ . 1 > , 1 ( n) = J I max < I 


where J(m ) is a Jordan bound, which is not greater than (\/8 m + l) ” m — (\/8 m — l)' 
Then 

d = ( K3 yw-\ ( 3 ) 


,2m 2 


It is well-known in the theory of linear algebraic groups that any algebraic subgroup of a 
diagonalizable group D is the intersection of kernels of some characters of D (see Proposition 
in the page 103 of Given a connected algebraic group H, the following proposition 

describes algebraic subgroups that are the intersections of some characters of H. 


Proposition 2.5 Suppose that H is a connected algebraic subgroup of GL n (Q). Then G is 
the intersection of kernels of some characters of H if and only if H is a proto-group of G. 


Proof. Assume that H is a proto-group of G. Let Xh • • • ,Xe be generators of X(H), the 
group of characters of H. Define a map if : H — >• (Q X ) £ given by if(h) = (xi (h), ■ ■ ■ ,Xi(h)), 
where Q x denotes the multiplicative group of Q. Then if(H) is a diagonalizable group and 
if(G) is one of its algebraic subgroups. Due to Proposition in the page 103 of pi)], ip(G) is the 
intersection of kernels of some characters of if(H). Denote these characters by x\, ■ ■ ■ ,xi. 
Notice that if induces a group homomorphism 

if*:X((® X Y)^X(H) 

X^x°f>- 

We claim that G = C l i=1 ker (if*(xi)). Obviously, G C C\ l i=l ker (if*(xi)). Suppose that 
h € n( =1 ker (if*(xi)). Then Xi{if{h)) = 1 for all 1 < i < l. This implies that if(h) € if(G). 
Lemma B.9 of m states that H u = ker^). Hence ker(-0) C G and then h £ G. 

Conversely, G is the intersection of some characters of H . Then H u = ker^) C G. Since 
H u is connected, H u C G°. Thus H is a proto-group of G. □ 


The connection between proto-groups and u-ideals in k[Y, 1/det(T)] is the geometric objects 
so called /c-torsors, which are introduced in the next section. 
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2.2 fc-Torsors 


We shall use Zero(J) to denote the set of zeroes of J in GL n (k), where J is a subset of 
k[Y, 1/det(y)]. Suppose that Z C GL n (k) is a variety defined over k. We shall use 4(Z) to 
denote the set of all polynomials in k[Y, 1/det(Y)] that vanish on Z. 

Definition 2.6 (see Definition 3.13 of f£2j/ ) Let Z C GL n (k) be a variety defined over k 
and H an algebraic subgroup of GL n (k) defined over k. Z is said to be a k-torsor for H if 
for any z\,z^ G Z, there is a unique h € H such that zi = Z 2 h. A k-torsor Z for H is said 
to be trivial if Z n GL n (k) 0, i.e. Z = BH for some B € GL n {k). 

Let / be a maximal c-ideal of k[Y, 1/det(Y)]. Then one has that 

Proposition 2.7 (Proposition 1.20, page 15 of \23j) Zero(I) is a trivial k-torsor for stab(I). 

Suppose that H is a connected algebraic subgroup of GL n (k), which is defined over Q, and 
Z is a trivial k -torsor for H. Then for any B G Z G GL n (k), the map given by 

k[Y, l/det(y)]/4(if) -> k[Y, 1/det(Y)\/I k (Z) (4) 

p(y) -> p(b~ 1 y) 

is an isomorphism of fc-algebras. A theorem of Rosenlicht ( [1411181 [21 j) implies that invertible 
regular functions on Z are closely related to characters of H. This theorem states: let H 
be a connected linear algebraic group defined over an algebraically closed field k and y be a 
regular function on H with 1/y also a regular function, then y is a k multiple of a character. 
Notice that characters of H can be viewed as elements in Q[T, 1/det (Y)]/Iq(H). 

Lemma 2.8 Suppose that J is a prime a-ideal of k[Y, 1/det(y)] and Zero(J) is a trivial k- 
torsor for H. Let B G Zero(J)nGL n (fc). Ifx is a character of H, then x(B~ l Y) is invertible 
in k[Y, l/det(y)]/J. Conversely, if P is an invertible element in k[Y, 1/det(y)]/ J, then 
P = rx(B~ 1 Y) for some r G k and some character x of H. 

Proof. We only need to prove the second assertion. Since Q is algebraically closed, H 
viewed as a linear algebraic group defined over k is still connected. The map <[4|) implies that 
P{BY) is invertible in k[Y, 1/ det(y)]/4(iL). Applying the above theorem of Rosenlicht to 
P(BY), one has that P(BY ) = ry for some r G k and some character y. Observe that 
k[Y, 1/det(y )]/J is a e-extension ring of k. Due to Lemma 1.19 in the page 15 of [23] . 
( k[Y, 1/det(y)]/J) n k = k. Hence (k[Y, 1/det(y)]/4(iL)) 0 k = k. We then conclude that 
r G k and P = rx(B~ l Y). □ 

In Section 4, one will see that invertible elements of k[Y, 1/det (Y)\/J are actually a- 
hypergeometric over k. In the case that J is a proto-maximal u-ideal, algebraic relations 
among these cr-hypergeometric elements will reveal the characters of H that determine the 
Galois group G. 
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2.3 Some related problems 

In this paper, we shall need the algorithmic solutions of the following problems. 

(PI) Given an ideal in k[Y], compute a Grobner basis of it with respect to some monomial 
ordering. The reader is referred to Section 2.7 of [1] and Section 5.5 of [Tj for the 
algorithms. 

(P2) Given an unmixed ideal in k[Y], compute its radical and its associated primes. There 
are several methods for this problem, for instance the methods presented in [8], Section 
4 of [6| . Section 8.7 of |JLJ, parts 36 and 42 of m- 

(P3) Compute the Galois group of linear difference equations of diagonal form. Equivalently, 
given ,be € k, compute a set of generators of the following Z-module: 

l 

3/e/cs.t. n h T = 

i =1 

When k = Q(x), a method was described in Section 2.2 of |23| . Using the results in 
Section 3.2 of [5], one can adapt the method in (23] to solve the problem with k = Q(x). 
This problem is the bottleneck in extending our algorithm to the equations over a larger 
basic held. 

(P4) Give linear difference equations with coefficients in k , compute all hypergeometric 
solutions. The reader is referred to ([HI [16]) for algorithms. 

3 The computation of a proto-maximal tr-ideal 

Let P be a fundamental matrix of (P) and let d be a positive integer or oo. Denote 

lF,d = {{P(Y) € k[Y]< d | P(P) = 0}), (5) 

where k\Y]< d denotes the set of polynomials in k [Y ] with degrees not greater than d, and 
(*) denotes the ideal in k[Y, 1/det(T)] generated by *. When d = oo, Ip td is equal to If 
that is defined in Introduction. One can readily verify that Ip d is a c-ideal and furthermore 
Ip is a maximal cr-ideal. The fact that k[Y, 1/det(T)] is a noetherian ring implies that for 
sufficiently large d, Ip d is a proto-maximal a-ideal. Therefore to achieve a proto-maximal 
cr-ideal, one only needs to solve the following two problems: (a) Given an integer d, how to 
compute If^I ( b ) When is the integer d large enough such that Ip d is proto-maximal? 

3.1 The computation of Ip )d 

In [12], Kauers and Zimmerman presented an algorithm for computing generators for the ideal 
of algebraic relations among solutions of linear difference equations with constant coefficients. 
Their algorithm relies on the fact that one can explicitly write down solutions of the equations 
of such type. Here, our task is different. We only compute the ideal generated by algebraic 
relations with bounded degree, while we are interested in linear difference equations with 
coefficients in k. 


i z i) € ^ 
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We first show that which fundamental matrix F we take in this section. Let Sq be the 
difference ring of germs at infinity of Q (see Example 1.3 in the page 4 of [23] for the definiton). 
Let p be a nonnegative integer such that i is not a pole of entries of A and det(A(i)) /Off 
i > p and Z p € GL n (Q). Define an element of GL n (<SQ), say Z = (Z 0 , Z 1: ■ ■ ■), as follows: 
Zi = 0 for 0 < i < p — 1 and Z l+ \ = A(i)Zi for i > p. Define a map 

/ : k[Y, 1/ det(y)] -b Sq 


as follows: 

for / € k, 'ip(f) = (0, • • • ,0, f(i), f(i + 1), • • • ,) and = Z 

where i is a nonnegative integer such that j is not a pole of / if j > i. Proposition 4.1 in 
the page 45 of [23] states that ip induces an embedding of k[Y, 1/det(Y)]// into <Sq, where 
I = ker(/) that is a maximal cr-ideal. Let F be the image of Y in k[Y, 1/det(Y)]//. From 
this construction, we have that I fa = Iz,d- 

The results in Appendix A imply that one can compute an integer l such that Ipd has 
a set of generators consisting of polynomials in Q[x][Y] whose degrees in x are not greater 
than l. Let N = ) — 1 and mo, • • • , mjy be all elements in Z"q with |mj| < d. Write 

P = ('j(t.+ 1 )+i Y mj for P € If,cI, where Y mi = WyJ’l' 0,1 with = (m t jj). We can then 

reduce the original problem to the following problem: find a basis of the vector space 


-1 

1 (co,ci, • • 

l N 'j 

•,c (JV+1)(WH )GQ (iV+1)(<+1) 


l 

ir 

o 

II 

o 


We are going to solve the latter problem. Observe that a ) is a /s-linear combination 
of the monomials F m °, • • • , x l F m i , • • • , x^F mN . Hence there is a nonzero linear difference 
operator L in Q[x][<9] such that L(x i F m - 7 ) = 0 for all i with 0 < i < £ and 0 < j < N. 
This operator L can be computed using the equation ([1]). Notice that at present, we do not 
know the ideal I and thus do not know F. Fortunately, one can easily compute the sequence 
solution Z, which can be considered as a difference analogue of formal power series solutions 
of linear differential equations,. 

For the convenience, write (H]) and L in the form of linear recurrence equations 


Y m+ 1 = A(m)Y m , m > p 


( 6 ) 


and 


L = ai(m)y m+ i + a;_i(m) 2 / m+i _i 4-b a 0 (m)y m , m>v 


(7) 


where p is a positive integer such that i is not a pole of entries of A(x) and det(A(i)) / 0 
for all i > p, and u is an integer greater than integer roots of ai(x)ao(x) = 0. One easily sees 
that 


Lemma 3.1 Assume that {s u ,s u+ 1 ,--- ,} is a solution of If there is a nonnegative 
integer j such that s u+ j = • • • = s u+ i-\+j = 0, then Si = 0 for all i > u. 
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Let k be an integer greater than p and v. Notice that the sequence {Z p , Z p+ i, ■ ■ • ,} is a 
solution of dU) and for all 0 <i<£ and 0 < j < N, the sequence , (« + 1)®Z^ 1( ■ ■ ■ , } 

is a solution of (0. Set 

e n 

c j(e+i)+i xt Y mj , w here c — (co,-- - , C(jv+i)(£fi)-i) e ^ ^ 

i =0 j =0 

Then the sequence {P c (k, Z K ), (P c (k + 1, Z re+ i), • • • , } is also a solution of 0. 

Proposition 3.2 c E U if and only if P c (i, Zf) = 0 for all k < i < k + l — 1. 

Proof. Assume that c G U. Then P c (x,F) = 0 and thus i/j(P c (x, F)) = 0. In other words, 
there is a positive integer j such that P c (i,Zi) = 0 for all i > j. Lemma 13.11 implies that 
P c (i,Zi ) = 0 for all n < i < k + l — 1. Conversely, suppose that P c (i,Zi) = 0 for all 
k < i < k + 1 — 1. By Lemma 13.11 again. P c {i,Zf) = 0 for all i > k. This implies that 
if(P c (x, F)) = 0. Equivalently, P c (x,F) = 0. Hence c G U. □ 


The conditions P c (i, Zf) = 0 for all«;<i<K + / — 1 induce a linear system for c. Solving 
this system, we obtain a basis of U. 

Algorithm 3.3 Compute a basis of 

(i) Using the results in Appendix A, compute an integer t such that Ip ^ has generators 
consisting of polynomials in Q[x][T] whose degrees in x are not greater than i. 

(ii) Construct a nonzero operator L in Q[x][<9] that annihilates sT mj for all 0 < * < i and 
0<j<N, where mo, • • ■ , mjy are all elements in Z> 0 satisfying |mj| < d 

(iii) Let k be an integer that is greater than both p and all integer roots of the leading and 
trailing coefficients of L. 


(iv) Compute Z K , Z K+ 1, ■ ■ ■ , Z K+ i_i, where l = ord(L). Set 

t N 

P c (x,Y) = EE Cj(i+l)-\-iX Y J , C ( Cq , • • • , ) . 

i=0 j=0 

Putting P c (k , Z K ) = • • • = P c (k + 1 — 1, Z K+ i_\) = 0, we obtain a linear system L in 
co, ci, ■ ■ ■ , C(jv+l)(m)—i- 

(AT+1)(£+1) 


}■ 


(v) Solve C and return <^Pc(x,Y) c € Zero(£) n 

Example 3.4 Consider the Fibonacci numbers F[n). It satisfies that 

F(n + 1)\ _ /0 l\ / F(n) 

F(n + 2) i ll 1 J \F(n + 1) 


Z = 





Let 



We are going to calculate Iz. 2 - Using the results in Appendix A, one sees that there are gen¬ 
erators of Iz ,2 whose degrees in x are zero. Let mo, • • • ,mi 4 he all vectors in Z 4 0 satisfying 
|mj| < 2. Let 

L = d 6 — Ad 5 + 2 d 4 + 6 < 9 3 - Ad 2 -29 + 1 . 

Then L annihilates Z m * for all 0 < i < 14 . Computing the first 6 terms ofZ , denoted by Z, 
for i = 0, • • • , 5 . Set c = (co, ci, ■ ■ ■ , C14) and let P c {x, Y) he defined as in the step ( d ). Then 

P c (x, Z 0 ) = c 0 + Cl + c 4 + c 5 + c 8 + C14, 

P c {x, Z\) = Co + C2 + C3 + C4 + Cg + C10 + Cu + C12 + C13 + C14, 

P c (x, Z2) = Co + Cl + C2 + C3 + 2C4 + C5 + C6 + C7 + 2 c 8 + Cg + C10 + 2 cn + C12 + 2ci3 + 4ci4, 
Pc(x , Z3) = Co + Cl + 2C2 + 2C3 + 3 c 4 + • • • + 4ci2 + 6ci3 + 9 C 14 , 

Pc{x , Z/f) = Co + 2 ci + 3 c 2 + 3 c s + 5 c 4 + ■ • • + 9ci2 + 15 ci 3 + 25 ci 4 , 

P c (x, Z 5 ) = Co + 3 ci + 5 c 2 + 5 c 3 + 8C4 4 -f 25 ci 2 + 40 ci 3 + 64 ci 4 . 

Solving the linear system {P c (x, Zf)\i = 0 , ■ ■ ■ , 5 }, one has that 

Co = 0 , Cl = —C4, C2 = —C4 — C3, C5 = —C 8 — C14, 

C6 = —c 8 — 2 cm — C7 — C11 — C13, C9 = —C14 — C10 — C\\ — C\2 ~ C13. 

From this, one sees that /z,2 is generated by 

2/2,1 — 2/1,2, 2/2,2 — 2/1,2 — 2/1,1- 

3.2 When is Ip,d proto-maximal? 

Let d be as in ([ 3 ]). In this section, we shall show that I F j is proto-maximal. Before proving 
this, we first describe some properties of I F d- Note that Ipd is contained in a maximal 
cr-ideal /. Proposition 1.20 in the page 15 of ([ 23 ]) states that Zero(/) is a trivial /c-torsor for 
stab(J). We show that a similar property holds for Ip^, i.e. Zero (lF,d) is a trivial /s-torsor 
for stab(I_F,d)- As Zero(I) is a trivial fc-torsor for stab(J), Zero(I) n GL n (fe) / 0 . Therefore 
Zero (lF,d) H GL n (A:) / 0 . For short, we denote by Hp,d the stabilizer of If./I- 

Lemma 3.5 Let B he an element of Zevo(Ip^) G GY n (k). Then 

I F ,d = ({QiB^Y) | Q G I^(H F4 ) nQ[Y]<d}) ■ 

Proof. Denote by J F the right-hand side. Suppose that P is an element of k\Y]<d with 
P(F) = 0 . Then for each h G P(Yh) € I fa and therefore P(Bh) = 0 . Write 

i 

P{BY) = Y J C i P i {Y) 

2—1 

where Pi(Y) G Q[P] and ci, ■ ■ ■ , c; are linearly independent over Q. Obviously, for all i with 
1 <i < l, the degree of Pi(Y) is not greater than d and Pi(h) = 0 for all h G H F d- In other 
words, Pi(Y) G Iq(H F4 ) 0 Q[T]<d for all i = 1 , • • • , 1 . Hence P G Jb and then I fa Q Jb- 
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Notice that Ip is a maximal er-ideal that contains I fa- Let G = stab(I^). Observe that 
the action of GL n (Q) on k[Y] preserves the degrees of polynomials. From the definition 
of If,<1i one sees that G C Hfa- Due to Proposition 12.71 Zero (Ip) = BG(k ) for any B G 
Zero (Ip) 0 GL n (fc). This implies that Zero (Ip) C Zero(J^) and thus Jg C Ip. As F is 
a zero of Ip, it is also a zero of Jg. This together with the fact that Jg is generated by 
polynomials in k[Y]< d implies that Jg C Ip d . On the other hand, the previous result shows 
that If a C Jg. Consequently, I fa = Jg- It remains to prove that Jp = Jg. 

We first have that Jg C Jp. Define a ^-automorphism f of k[Y] as follows: 

f{P{Y)) = P(BB~ l Y). 

Then f(Jg) = Jp , which implies that Jp C f(Jp) C f 2 {Jp ) C ■ ■ ■. Due to the noetherian 
property of k[Y, 1/det(Y)], Jp = f(Jp). In the sequel, Jp = Jg. □ 

The above lemma has the following corollaries. 

Corollary 3.6 Let B be an element ofZero(Ip t d) O GL n (fc). Then 

Zero (I Ftd ) = BHp td (k) 

i.e. Zero (If^) is a trivial k-torsor for Hp t d. 

Corollary 3.7 Let I uv be an associated prime of Ip d . Then stab(/j rr ) = Hf d . Moreover 
Zero(/i rr ) is a trivial k-torsor for H° Fd . 

Proof. Let B be an element of Zero(/i rr ) n GL n (fc). By Corollary 13.61 

Zero(/i rr ) n GL n (k) = BHifk) 

where Hi is an irreducible component of Hp d - Since B G Zero(/i rr ), Hi = Hf d . Thus 
Zero(/i rr ) is a trivial fc-torsor for H Fd . For each g G GL n (fe), one can define an isomorphism 
4> g of GL n (L) given by </> g (Z) = Zg. As Ii„ is prime, one can verify that h G stab(/i rr ) if and 
only if <fh(BHp d (k)) = BH Fd (k). On the other hand, for h G GL n (Q), (fh(BH° Fd (k)) = 
BH Fd (k) if and only if h G H° Fd . Therefore stab(/ irr ) = H° Fd . □ 

Lemma 3.8 Suppose that F is a fundamental matrix of W- If F is a zero of Ip d , then 

I F,d = k F d- 

Proof. From the assumption, one has that Ip d C I F d . Observe that F = Fh for some 
h G GL ri (Q). The definition of Ip )d implies that < f>h{^Fd ) = IpcF w Lere the homomorphism 
fh is given by fh(Y) = Yh. The successive applications of fh to If^ L=. Ip d induce that 

Ip )d C fh(Ip,d ) ^ ^ 4>h(lF,d ) ^ ‘ ‘ ‘ • 

The noetherian property of the ring k[Y, 1/det(Y)] implies that I F d = fh(Ip d)- Therefore 
fh(lF,d ) = fh{Ip,d) > i- e - IfA = Ip,d- Dl 
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Proposition 3.9 stab(/p d is a proto-group of stab(J), where I is any maximal a-ideal 
containing I F g- In other words, I F i is proto-maximal. 

Proof. Let G = stab(I) and H be an algebraic subgroup of GL n (Q) that is bounded by 
d and is a proto-group of G. Such H exists by Proposition 12.41 Observe that there is a 
fundamental matrix F such that I = Ip. Since I F p C /, we have that Ip r = I F j due to 
Lemma 13.81 Let B be an element of Zero(J) n GL„(fc) and let 

J = { Q(B~ l Y ) I Q € %(P) flQ[y]<j} . 

Since H is bounded by d, Zero (J) = BH(k). By Proposition 12.71 Zero(I) = BG(k). There¬ 
fore J C /, because I is radical and H is a proto-group of G. Note that F is a zero of I. 
One then has that F is a zero of J. This implies that 


J — IF,d ~ IF,d — I' 


The first inclusion holds because J is generated by a set of polynomials in A;[Y’] < r. Lemma f3.5l 
then implies that 

g<h f ^<h. 

Then the proposition follows from Remark 12.21 □ 


Example 3.10 Consider 


Using the method developed in Section 3.1, we can compute a a-ideal 



( 8 ) 


Iz,2 = (2/1,12/1,2) 2/1,12/1,3) 2/1,12/2,1, 2/1,12/2,3, 2/1,12/3,1, 2/1,12/3,2, 2/1,22/1,3, 2/1,22/2,1, 2/1,22/2,2, 2/1,22/3,2, 2/1,22/3,3 
2/1,32/2,2,2/1,32/2,3, 2/1,32/3,1, 2/1,32/3,3, 2/2,12/2,2,2/2,12/2,3, 2/2,12/3,1, 2/2,12/3,3, 2/2,22/2,3, 2/2,22/3,1, 2/2,22/3,2 
2/2,32/3,2,2/2,32/3,3, 2/3,12/3,2,2/3,12/3,3, 2/3,22/3,3)- 


Furthermore, one has that 


stab(/ Zi2 ) 


(a 0 0\ 

0^0 

\0 0 1 ) 


ah ± 0 > (J 


/ 0 a 0\ 

0 0 p \ 
\7 0 0 j 


ah ± 0 > (J 


(0 

0 


1 

[p 

0 

° 

ah ^ 0 

\o 

7 

0 / 

J 


Since all elements in stab(/z. 2 ) is semi-simple, stab(/z, 2 ) is a proto-Galois group of (1\) 
over k, i.e. Iz ,2 is a proto-maximal a-ideal. 


4 The computation of a maximal a '’-ideal 

The results in the previous section enable us to calculate a proto-maximal a-ideal. Suppose 
that we have obtained a proto-maximal a-ideal, say I F c j. Let /; rr be an associated prime of 
I F p It can be obtained by the algorithmic solutions of the problem (PI). Since I pFL is a 
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a-ideal, I nv is a a^-ideal for some positive integer 8. In the following, we will enlarge I ur to 
obtain a maximal a d -ideal. By Corollary 13.71 one sees that for any B G Zero(/j rr ) n GL n (7), 

Zero(/i rr ) = BH° F ^(k) and stab(Ii rr ) = 77° ~ F (9) 

Let Is be a maximal a 5 -ideal that contains 7i rr . Proposition 12.71 implies that Zero(Is) is a 
trivial 7-torsor for stab(/< 5 ), i.e. Zero (Is) = BG$(k ) where B G Zero (Is) 0 GL n (7). Then 
the equality @ induces that stab(Is) C H° F We shall show that H° F ~ is a proto-group of 
stab (Is). 

Lemma 4.1 Let I be a maximal cr s -ideal and I = I D a(7) D • • • n a s ~ l (I). Then 

(a) I is a maximal a-ideal. 

( b ) [stab(I) : stab(7)] < 8. 

Proof, (a). Suppose that I is a maximal a-ideal. Let J be a maximal a^-ideal containing 
I. Then I C n^ZQa*(J). On the other hand, nfz^a l (J) is a a-ideal and thus it is equal 
to I. For any g G GL n (Q), one can define a a-isomorphism <j)g of k[Y, 1/det(T)] given by 
4> g (Y ) = Yg. From the uniqueness of the Picard Vessiot extensions, one can easily see that 
there is g € GL n (Q) such that 4> g (J) = I. This implies that (f> g {I) = /. Hence I is a maximal 
a-ideal. 

( b). Let G = stab(I) and G = stab(I). Let B be an element of Zero(I) fl GL n (k). Due 
to Proposition 12.71 one has that 

Zero(I) = BG(k ) and Zero(I) = BG(k). (10) 

Meanwhile, all a*(/) are maximal a^-ideals. Hence there are gi,--- ,gs~ i € GL ra (Q) such 
that (pg i (a' l (I )) = I. This implies that 

Zero(a*(/)) = BG(k)gi , i = 0,1, • • • , 5 — 1. (11) 

The equalities (flOl) and (fTll) imply that G = U^Zq Gg\. In the sequel, [G : G\ < 8. □ 

Let I = Is O u(Is) O ••• O a* 5 ^ 1 ^). Then I Ff j C I. The above lemma together with 
Proposition 13.91 induces that H F ^ is a proto-group of stab(I), i.e. 

( H Fd) ^ (stab (J))° < stab (I) < H FF 

Observe that [H f B = \H° -) . Due to the above lemma again, (stab(7))° = (stab (Is))°■ 
V ,a / u V u 

Thus 

(■ H °F,d) u - ( stab ( J 5))° < Stab (Is) < H° Fd 

i.e. H° f ~ is a proto-group of stab (Is). Proposition [23] implies that stab(Is) is the intersection 

of kernels of some characters of H° This will enable us to construct Is- Suppose that 

F,d 

Xi, ■ ■ ■ ,Xi are characters of H° - satisfying 

r ,CL 

ker(xi) n • • • n ker(x^) = stab( 7 , 5 ). 

Then we have the following lemma. 
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Lemma 4.2 Let B be an element of Zevo(Is) n GL n (k) and 


S = I irr U{xi(B~ 1 Y)-l \i = !,-•• ,/}. 


Then Zero(Is) = Zero(S). 

Proof. Let G$ = stab(i,s). It suffices to show that Zero(S) = BGs(k). Observe that B G 
Zero(/i rr ) O GL n (/c), which implies that Zero(/i rr ) = BH° F ~(k). Suppose that Z G BGs(k). 
As Gs is the intersection of kernels of the characters X\, - ■ ■ , Xz, one sees that Z G Zero(S). 
Conversely, assume that Z G Zero(S). Then Z G Zero(/i rr ) and thus Z = Bh for some 
h £ H° F ~{k). Meanwhile for each i = 1, • • • , l, 


Xi{B 1 Z) = Xi(h ) = 1 . 

This implies that h G Gs(k). Therefore Zero(S) = BGs(k). □ 

Proposition 12.71 states that invertible elements of k[Y, l/det(T)]/ii rr are k multiples of 
characters of H° F ~. Precisely, let P be an invertible element of k[Y,l/ det(Y)\/I[ rr , then 

P = rx(B~ 1 Y) for some r G k and some character x of H° f By the above lemma, to com¬ 
pute Is, it suffices to find invertible elements of k[Y, 1/det(y)]//i rr that take constant values 
on Zero (Is). In the following, we first show that invertible elements of k[Y, 1/det(y)]/I; rr 
are actually o^-hypergeometric over k and then prove that algebraic relations among cr s - 
hypergeometric elements take constant values on Zero (I 5 ) and enable us to find Is- We start 
with a definition. 

Definition 4.3 A nonzero element P of k[Y', 1/ det(y)]//i rr is said to be -hypergeometric 
over k if P is invertible in k[Y, 1/ det(T)]/ I ar and a s (P) = rP for some r G k. 

Let P±, P 2 be two cr 5 -hypergeometric elements over k of k[Y, 1/det(!")]//;„. We say P\ and 
P 2 are similar if there is r G k such that P\ = rP 2 - 

Proposition 4.4 Let B be an element of Zero(/; rr ) n GL n {k) and x a character of H° Fd 
that is represented by an element of Q[Y, 1/det(y)]. Then x{B~ l Y) is a a 5 -hypergeometric 
element over k of k[Y, 1/ det(Y)]/L TT . Furthermore, if x 1 cmd X2 are two distinct characters, 
then xi{B~ l Y) and X2 (B~ 1 Y) are not similar. 

Proof. Obviously, x(^~ 1 Y) is invertible in k[Y, 1/det(T)]//; rr . We first claim that 

a s (B~ 1 )AsB G H° F jfk). 

For any Q G -), by (|U|), Q(B~ 1 Y) G I nr . As I ur is a cr^-ideal, one has that 

Q(cr s (B~ 1 )AsY) G /; rr - Since B G Zero (Ii rr ), Q(cr s (B~ 1 )AsB) = 0 . This proves the claim. 
Now for any h G H° F ~(k), 

xia^B-^AsBh) - X (cT 5 (B- 1 )A 5 B)x(B~ 1 Bh) = 0 . 
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This implies that 


xia^B-^AsY) - X (a s (B~ 1 )A S B) X (B~ 1 Y) e I m . 

In other words, 

a^xiB-'Y)) - x(ct 5 (B- 1 )A s B) X (B- 1 Y) € I irr , 

i.e. x{B~ l Y) is a o^-hypergeometric element over k of k[Y. 1/det(Y)]/Ii rr . This proves the 
hrst assertion. 

Now assume that xi(B~ 1 Y) — rx 2 (B~ 1 Y) e Ii rr for some r e k. Then for any h e 7/° ~, 

Xi(M = yi(B _1 £fr) = rx2(B~ 1 Bh) = rx2(h). 

Particularly, putting h = I n , one then has that r = 1. Thus yi = X 2 , a contradiction. □ 


Let H 2 be as in d2|). Proposition B.17 of [7j states that X(H^~) has generators that are 
represented by polynomials in Q[Y]< K2 . Denote 

P is <7''-hypergeometric over k in k\Y, 1/det(Y)l//i rr , 

P- rQ $ 7 irr , V r g k, V Q € TL \ {P} 

and 

P € X(Hy, 

P-Qth{H° F4 ), VQeX\{P) 

Then X is a set of generators of X(i7° ~). Furthermore, we have that 

Corollary 4.5 There is a bijective map between TL and X. 

Proof. Let B G Zero(/i rr ) n GL n (k). We define a map r from TL to X(H° 3) as follows: 
t(P) = x where x £ satisfies P — rx{B~ l Y ) € /; rr for some r € k. By Proposi¬ 

tion E3 for each P e "H, there is a character y such that t(P) = y, and such character is 
unique by Proposition 14.41 Hence r is well-defined. From the definition of TL, one sees that 
t is injective. We shall prove that t(TL) = X. As the map defined in Q is an isomorphism, 
one has that 

P(BY)/r- X (Y)eI k (H° F>d y 

Hence y(Y) can be chosen to be a polynomial in Q[Y]< K2 . That is, y e X. Therefore 
t{TL) C X. Finally, Proposition 14.41 implies that t{TL) = X. □ 



x = {Pe Q[Y ]< K2 


n = \ p e k[Y ]< K2 


Algorithm IB. II in Appendix B enables us to compute TL. Suppose that TL = {Pi, ■ ■ ■ , P u }. 
Let bj be the certifications of Pj, i.e. cr s {Pj) — bjPj £ T rr for all 1 < j < n. Set 


Z = { (mi, ■ 




£ IT 


3 fee 


s.t. 

3 =1 


°\f) 
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Z is a finitely generated Z-module. The solution of the problem (P4) allows us to compute 
a set of generators of Z. Assume that mi, • • ■ , are generators of Z and further suppose 


that 





IP? 




j =1 

where 

fi € 

k and m,; = (m^i, • • • , 


where 

m ^ 

LLL i > 

m~ are in Z> 0 and 

( m D T 

pm _ 

n;= 

1 Pj lj where m = (mi, • • 

• j m v ). 


<r S (fi) 

fi 


, n, write m* = — m 8 , 

the vector (Pi, • • • , P^) and 


V = 




* = !,■■■ , 1 * 


}) 


It is easy to verify that V is a cr d -ideal. Let Is be a maximal cr^-ideal containing V. Then 
Proposition 4.6 Zero(P) = Zero (Is), he. Is = VV. 


Proof. Let B be an element of Zero(Is) n GL n (k) and G$ = stab (/a). Then due to Propo¬ 
sition 12.71 

Zero(I^) = BG$(k). 

The discussion after Lemma l4~T1 states that H° - is a proto-group of Gs- By Proposition 12.51 

r ,a 

Gs is the intersection of kernels of some characters of H° Let A be the set of these 

F,d 

characters. Observe that A is a set of generators of X(H° F -). Suppose that \ G A. Then 


X = Y[r(P,r\ ( 12 ) 

i =1 

where ati € Z and r is defined as in Corollary 14.51 By Corollary 14.51 for each i = 1, • • • , 1/, 
there is r, € k such that 

r(P i )(P- 1 P) - ri Pi € / irr . (13) 

Lemma [4~2l implies that — 1 € Is- Denote by Y the image of Y in k [Y. 1/ det(P)]// 5 . 

Then x(P _1 F) — 1 = 0. This together with (fl2l) and (fl3l) induces that 

V 

nW(y)-i = o. ( 14 ) 


Applying a s to (ThO) . one has that 


i =1 


Combining (1141) and (115|) . one has that 




n 

i=l 


v S ( —OLi\ 

° Vi ) 

„-oti 


( 15 ) 
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Set a = (oq,--- ,a u ) € IT. Then a G Z. So there are integers z i,--- , z fl such that 
a = 2imi H-b m M . 

Let Z be an element of Zero('P). Then one has that P mi (Z) = /, for all 1 < i < n, 
because P m * (Z) / 0. By (fT2l) and (fl3l) again, 

V V 

x(B~ 1 Z) - 1 = HriP^iB-'Z) - 1 = P Q (Z) JJrf - 1 
1=1 1=1 

[1 v [l V 

=n p* >m> (z> n <■? - 1 =n ff n ■c - »• 

i= 1 i=1 i= 1 i=l 


This implies that the polynomial x(B 1 X) — 1 takes a constant value on Zero(P). Particu¬ 
larly, putting Z = B, one has that xiB^B) — 1 = TT-Li /, Zi 715=1 r f' — 1 = 0. In the sequel, 
xiB-'Z) -1 = 0 for all Z € Zero(P). Therefore 

Zero(T’) C Zero(I; rr U {x(B _1 T) — 1 | x €E A}). 

The former set contains Zero(I^) and the latter one is equal to Zero (Is) by Lemma 14.21 
Consequently, Zero('P) = Zero(I^). □ 

Suppose that V has been calculated. One can then compute \fV by the methods developed 
in (|6], Section 8.7 of p]) and / = \fV O a(\fV) 0 ■ ■ ■ D cr' 5_1 (v / P) by the algorithm presented 
in (Section 6.3, page 260 of [I]). Then the ideal I is a maximal cr-ideal by Lemma 14.11 

Example 4.7 (Example I o'. 1 0 \ continued) We have the following irreducible decomposition: 


lz ,2 


— (2/1,1) 2/1,2,2/2,2) 2/2,3) 2/3,1) 2/3,3) O (2/1,1,2/1,3) 2/2,1,2/2,2) 2/3,2) 2/3,3) 


n (2/1,2, 2/1,3) 2/2,1, 2/2,3, 2/3,1) 2/3,2) • 


Set /j rr = (2/i,i)2/1,2,2/2,2)2/2,3)2/3,1,2/3,3)- Then one can easily verify that I lvv is a a’ 1 ’-ideal and 


stab(/j rr ) 


/a 0 0 \ 

0/3 0 

\0 0 7 / 


a /?7 / 0 > . 


The group of characters of X(stab(/i rr )) is generated by 2 / 1 , 1 , 2 / 2 , 2 ) 2/3,3- Thus we only need 
to compute a 3, -hypergeometric elements in k[Y, 1/det(Y’)]/Ji rr which are represented by lin¬ 
ear polynomials in k[Y], By Algorithm, \B. 11 we hare that 2 / 1 , 3 , 2 / 2 , 1 , 2 / 3,2 are a 3 -hypermetric 
elements of k[Y, 1/det(y)]//j rr and further they are not similar in pair. Precisely, 

o’ 3 ( 2 / 1 , 3 ) = (a; + 2)yi, 3 , cr 3 ( 2 / 2 , 1 ) = xy 2 , 1 cr 3 ( 2 / 3 , 2 ) = (® + 1)2/3, 2 - 

An eas// calculation implies that the only element (mi, m 2 , m 3 ) in Z 3 such /ha/ 

x mi (x + l) m2 (x + 2) m3 = a 3 (f)/f 

/or some / € k is (0,0,0). This implies that T m is a maximal a 3 , -ideal. 
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5 The algorithm and an example 


We are now ready to present the algorithm for computing the Galois group stab (/), where 

I is a maximal er-ideal of k\Y, 1/det(P)]. 

Algorithm 5.1 Input: linear difference equations of the form 0). 

Output: the Galois group of 0) over k. 

(i) Compute a proto-maximal ideal I F j by Alaorithm \3.3[ 

(ii) Using algorithms for the problem (P 2), compute an associated prime of I F denoted 

by I[ rr . Compute a positive integer S such that I m is a a 5 -ideal. 

(iii) By Alaorithm \B. 1[ compute a s -hypergeometric elements in k[Y, 1/det(y)]//; rr that are 
represented by polynomials in k\Y]< K2 , and are not similar in pair. Denote them by 
Pi,-" , Pv- 

(iv) Let bi be the certificates of Pi, i.e. v s {Pi) — 6jPj € I nT where bi € k and i = 1, • • • , u. 
Using the method for the problem (P3), compute a set of generators of the following 
h-module 


Z = 


Oi, • • • ,z v ) £ IP 


3fek x ,s.t. n b? 

i =0 



Denote those generators by mi, • • • ,m M . 

(v) Set P = (Pi, • • • ,P V ) and find f),, the element in k satisfying P m ' = & S (fi)/fi where 
i = 1, • • • , v. Set 


V = I irr U { P m » + - /. t P r 


* = !,••• ,/i 


}• 


where m^m^ are elements in Z> 0 satisfying mt — m 8 = m,; and m^ (m,- j 7 = 0. 

(vi) By the algorithms for the problem (PI) and the algorithm presented in (Section 6.3, 
page 260 ofW), compute Vv and 


i=vv rK^n-rv-H' 


(vii) Return stab(J). 

The correctness of the algorithm comes from the results presented in the previous sections. 


Remark 5.2 One may suspect that the complexity of the algorithm would be very high, 
since the integers d and k 2 given in m and (0) are quite large. These integers guarantee 
the terminate of the algorithm. However, as shown in Examvles \3. 1 dl and 4 . 7, these integers 
may be much larger than those required in practice. 


In the following, we give an example to illustrate the algorithm. 
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Example 5.3 Consider the following linear difference equations 

(vi\ (0 1 0 

a h/2 = \x o o 

W \° 0 h 

(i) Using the method developed in Section 3.1, we compute an ideal I generated by poly 

nomials in If ,2 n Q[E]: 

I = (2/3,2,2/3,1,2/2,3,2/2,12/2,2,2/1,3,2/1,22/2,2,2/1,12/2,1,2/1,12/1,2) • 

I is a a-ideal and 


( a 0 0 N 
stab(7) = ^ ( 0 /3 0 

,0 0 7y 


'0 a 0 N 


^ 0 f U 1 ( ^ 0 0 

0 0 7 , 


a /37 + 0 > • 


As all elements in stab(J) are semi-simple, stab(I) is a proto-maximal cr-ideal and thus 
I is a proto-maximal a-ideal. 

(ii) I is radical and one can compute its irreducible decomposition as follows: 

I = (2/1,1,2/1,3, 2/2,2, 2/2,3, 2/3,1, 2/3,2) n (2/1,2, 2/1,3, 2/2,1, 2/2,3,2/3,1, 2/3,2)• 

Set 7 irr = (2/1,1,2/1,3,2/2,2,2/2,3,2/3,1,2/3,2). Then 7 irr is a a 2 -ideal and 

stab(7 irr ) = {diag(a,/ 3 , 7 )|a /37 7 ^ 0}. 

(iii) Observe that the group of characters of stab(/i rr ) is generated by linear polynomials. 
Using Alaorithm \B.l[ we can find that a 2 -hypergeometric elements of k[Y, 1/det(y)]/; rr 
that are represented by linear polynomials in k[Y ] are 2/1,2,2/2,1,2/3,3- Precisely, 

1 


v (2/1,2) = ®2/l,2, cr (2/2,1) = (® + 1)2/2, 1, ^'(2/3,3) = 


(iv) 5e/ 


^ = < (zi,z 2 ,z 3 ) € Z 3 


3 f G k x ,s.t. x Zl (x + 1) 


22 


1 


x(x + 1) 


Z3 <?\f) 

f 


One sees that Z is generated by (1,1,1). 


(v) Let V = (/in- U {2/1,22/2,12/3,3 ~ 1})- One has that V is a radical ideal and thus is a 
maximal a 2 -ideal. 

(vi) Compute I = V n cx('P). One has that 

I = (2/3,2,2/3,1, 2/2,3, 2/2,22/2,1, 2/1,3,2/2,22/1,2, 2/1,22/1,12/3,3 — 2/2,1,2/1,22/2,12/3,3 — 2/1,2, 
2/1,22/2,12/3,3 + 2/1,12/2,22/3,3 — 1,2/1,12/2,1, 2/1,12/1,2)- 

(vii) Using the Grobner base computation, we have that 


'a 0 O' 
stab (7) = •( ( 0 (3 0 
0 0 7 ; 


0 a O' 


a /3 7 = 1 f U 1 | ^ 0 0 

0 0 7 j 


a/^7 = 1 
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A Coefficient bounds for generators of Ip^ 

Note that If 4 is generated by 

§ = { P (Y) € k[Y]< d | P(F) = 0}, 

which is a fe-vector space of finite dimension. We are going to find coefficient bounds for §. 
Precisely, we shall find an integer £ such that there is a basis of § satisfying that coefficients of 
elements in this basis are of degree < £. Let N = ( d+ d ) and M\, • • • , Mjv be the monomials 
in entries of F with degrees not greater than d. Observe that for a basis of §, it suffices to 
find a basis of the following vector space 


(ai, • • • , ajsr) €E k 


N 


N 


a%Mi = 0 


2=1 


Furthermore, one sees that (Mi, • • • , Mn) t is a solution of linear difference equations, which 
can be constructed from ([!]). Hence our original problem can be reduced to the following. 

Problem A.l Assume that v = (iq, • • • ,v n ) T is a nonzero solution of d, where the Vi are 
in some Picard-Vessiot extension ring ofk. Set 


W = {(a 1 ,a 2 , • • • , a n ) € k n \ a\V\ H-b a n v n = 0}. 


Find an integer £ depending on n and A, such that W has a basis consisting of vectors whose 
entries are of degree not greater than £. 

Without loss of generality, we may assume that v\, • • • ,v r are linearly independent over k 
and 

Vr-\-i — C^ ] Vi “b "P Ci^Vxi 2 — 1, , Cl r . 

For all i with 1 < i < n — r, denote c* = (cj ; i,--- , , Cj >n ) where Ci, r+ i = —1 and 

Ci )T+ j = 0 for 1 < j < n — r and j / %. Then {ci, • • • , c n _ r } is a basis of W. Actually, for 
any a = (ai, • • • , a n ) € W, we have that a = — (a r +iCi + • • • + a n c n _ r ). In the following, we 
are going to find a bound for deg(cjj), where i = 1, • • • , n — r, j = 1, • • • , r. Let V be the 
solution space of (P) and 

V = {w € V | c,;w 7 =0, V i = 1, ■ ■ ■ ,n — r}. 

Then V is a Q-vector space of finite dimension. Moreover, we have 

Lemma A.2 dim(H) = r. 

Proof. Clearly, V € V. Suppose that {vi, • • • , v^} is a basis of the vector space over Q 
spanned by the orbit of v under the action of Gal(K/k), the Galois group of (P), where K is 
the ring of fractions of the Picard Vessoit extension of k for P). Then Vj € V for all i with 
1 < z < //. In the sequel, dim(P) > /a. In the following, we shall prove that fj, > r. Denote 
the matrix consisting of the first /r rows of (vi, • • • , v^) by D and the remaining one by U. 

For any € Gal(A/A:), there is [<j>\ G GL M (Q) such that <j>(D) = D[<f\ and <j)(U) = U[(j)\. 
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Without loss of generality, we may assume that det (D) / 0. As for any 4> £ Gal (K/k), 
4>(det(D)) = det(D) det([^6]). One sees that det(D) is invertible in K and therefore D is 
invertible. Now for any 4> € Gal(A"/fc), 

${UD~ l ) = U&M-'D - 1 = UD~\ 

The Galois theory implies that C = UD~ l € Set C = (C, I n -^)- Then 



Particularly, CV = 0. This implies that dim(VP) = n — r > n — /i and then n > r. So 
dim(P) > r. On the other hand, one has that dirn(y) + n — r < n and then dim(P) < r. 
Hence dim(H) = r. □ 

Assume that {vi = v, V 2 , • • • , v r } is a basis of V and M is the n x r matrix consisting of 
the vectors vi, • • • , v r . For 1 < i\ < ■ ■ ■ < i r < n, denote the determinant of the sub-matrix 
consisting of the ii-th,Z 2 -th, ••• , z r -th rows of M by du,i 2 ,--- ,i r - Then an easy calculation 
implies that 

= ■■■ ,r ) where ,i r € k. 

In particular, 6i, 2 ,...,j-i,j+i,...,r,r+i = (-1 Let b = >r ,y , 6„_ r+ i )n _ r+2 ,...,n) T - 

On the other hand, one can construct from A an invertible matrix A r with entries in k such 
that hdi£,— ,r is a solution of cr(Y) = A r Y. Notice that the matrix A r only depends on A 
and r. Moreover, one can easily verify that di } 2 ,— , r is hypergeometric over k. This implies 
that bc?i j2) ... )T . is a hypergeometric solution. By means of cyclic vector, the system of the 
form dU) can be reduced into a scale linear difference equation. Then algorithms developed 
in (ia m) allow us to find all hypergeometric solutions of (JTJ). Therefore one can find an 
integer 1/2 such that hypergeometric solutions of cr(Y) = A r Y are of the form w h where 
h is hypergeometric over k and w is a vector whose entries are elements in k with degree 
not greater than i/2. Particularly, bc?i )2) ... )7 . = w h where w = (w i, • • • ,w n ) € k n satisfying 
deg(u)j) < i/2 and h is hypgeometric over k. Observe that = 1- Then one has that 

b = w/wi. Hence entries of b are of degree < £, i.e. deg(cjj) < i. 

In the case that we do not know the dimension of V, we can take r = 1, 2, ■ ■ ■ ,n and 
construct the corresponding systems <j(Y) = A{Y, ■ ■ ■ ,cr(Y) = A n Y respectively. Compute 
all hyper geometric solutions of these systems and let £/2 be an integer such that these 
hyper geometric solutions are of the form w h where h is hyper geometric over k and w is a 
vector whose entries are rational functions in x with degrees not greater than 1/2. Then we 
have that deg(cjj) < i. This solves Problem A.l. 

B tx 5 -Hypergeometric elements 

We shall describe a method to compute cr^-hypergeometric elements in k[Y, 1/det(y)]//i rr . 
In fact, we are not going to calculate all crMiypergeometric elements in k[Y, 1/det(y)]//i rr . 
Instead, we only find those (j 5 -hypergeometric elements that are represented by polynomials 
in k\Y] with degrees not greater than d and furthermore that are not similar in pair. Assume 
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that mi,-- - , are polynomials in k\Y]<d satisfying that {mi,-- - , m^} is a fc-basis of 
M*1<dMrrn k\Y]<d)i where rhj is the image of m,. By the Grobner base computation, one 
can find these rrij. As a s preserves the degrees of elements of k[Y], there is A G GL fik) such 
that 

^((mi,%-" ,*<)) = (mi,m 2 ,--- ,m e)A. 

The invertible matrix A can be constructed from A. Now suppose that P = Y) c^mj is a 
er^-hypergeometric element, where Cj G k, i.e. c s (P) — rP € I\ rr for some r G k. Then one 
can verify that ci, • • • , eg and r satisfying 


/ ci N 


Aa° 


'cA 


= r 


\C£/ \ClJ 

Let h be a cr 5 -hypergeometric element satisfying cr 5 (h) = rh. Then (ci,-- - ,cg) T h is a a s - 
hypergeometric solutions of the following linear difference equations 


a 5 (Y) = A~ 1 Y. 


(17) 


Consequently, for those ci, ■ ■ ■ , eg and r, it suffices to find all cr d -hypergeometric solutions of 
the above linear difference equations. The algorithms for computing all <r 5 -hypergeometric 
solutions of m can be found at m Il6]). Particularly, one can find (7 5 -hypergeometric 
solutions ci/ii, • • • ,cihi that are not similar in pair where h\, ■ ■ ■ ,hi are a 6 -hypergeometric 
and ci, • • • , ci are vectors with entries in k. Here two vectors hi, h 2 are said to be similar 
if hi = rh 2 for some r G k x . Furthermore, if h is a cr^-hypergeometric solution of (1171) . 
then there is a unique j with 1 < j < l satisfying h = bcjhj for some b G k. Write 
c i = (cj,i, • • • , Ci/) and set P { = Yf/mi where 7 = 1,2,--- , l. Then a 5 (Pi) - riPi G / irr 

for some r* G k. It remains to select those Pj that are invertible in k[Y, 1/det(y)]//; rr . Note 
that Pi is invertible in k[Y, 1/det(y)]//; rr if and only if Zero (Pi) n Zero(/i rr ) = 0. The latter 
condition can be detected by the Grobner base computation. Precisely, it suffices to decide 
if 1 is in the ideal (/i rr) Pi) ■ The previous results are summarized in the following algorithm. 


Algorithm B.l Compute all a s -hypergeometric elements in k[Y, 1/det(y)]//; rr that are 
represented by polynomials in k[Y]<d end are not similar in pair. 

(а) Compute a Grobner basis for I m D k[Y] and then find the monomials mi, ■ ■ ■ ,m^ in 
k[Y]<d such that {mi, • • • ,m£} is a k-basis of k[Y]<d/(fivrCk[Y])<d, where rh, denotes 
the image of m; in k[Y, 1/ det(T)]//;„. 

(б) Construct an invertible matrix A G GL i(k) such that 

^((mi, m 2 , • • • , m^)) = (mi, m 2 , • • • , m^)A 

(c) Compute a s -hypergeometric elements 

a s (Y) = A~ X Y 

that are not similar in pair, say C\h\, ■ ■ ■ , C( h/ , where h\, ■ ■ ■ ,hi are a 5 -hypergeometric 
and ci, • • • , ci are vectors with entries in k. 
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(d) Write c. t = (a t i, • • • , Ci/) and set Pi = Ylj =1 where i = 1, 2, • • • ,1. 

(e) Decide whether k[Y, z\ = (I- ur P\k\Y], Pi,det(Y)z — 1) by the Grobner base computation. 
Return those Pi satisfying (I- ur n k[Y], Pi,det(Y)z — 1) = k[Y, z]. 
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